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Abstract 

Hernandez and Smirnov discovered an interesting formula to parametrize each column of a 
neutrino mixing matrix by six integers related to the residual symmetry. We point out that these 
six integers are not independent, and propose a way to find the allowed combinations using structure 
constants of finite groups. 
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I. INTRODUCTION 



Symmetry of leptonic mixing can be studied from the top down, or from the bottom up. In 
the former case, a symmetry group Q is introduced from which the leptonic mixing patterns 
are computed. In the latter case, Q is derived from the residual symmetries compatible with 
experimental mixing. In this article we discuss the latter. 

It was recently found [1, 2] that each column of a neutrino mixing matrix can be 
parametrized by six integers related to residual symmetry. This is a great step forward 
because this formula for mixing vectors is explicit and easily calculable. In contrast, it is 
considerably harder and much more lengthy to produce mixing vectors from the top down, 
because in that case many eigenvectors of many group elements of many 3-dimensional 
irreducible representations of many groups have to be computed [3, 4]. 

Nevertheless, there are hidden difficulties with this bottom-up apprach. As we shall see in 
Sec. II, without a constraint between these six parameters, there may be no group containing 
the intended residual symmetry, much less a finite one with a moderate order. The purpose 
of this article is to find out the constraints between these parameters. 

After a review of the main formula [1] in Sec. II, and some of its consequences, we will 
embark on a way to search for the constraints in Sec. III. We need to do a scan to find the 
allowed parameters, but this scan is much simpler than the top-down scan [3, 4] because 
only class functions are needed here rather than the individual group elements. Examples of 
the search are given in Sec. IV, with the help of two GAP-language [7] programs discussed 
in Appendix A. 

II. REVIEW OF THE PARAMETRIZATION 

We follow the notations of [1]. The symmetry generator of Z m in the charged-lepton 
sector is denoted by T, rather than F as in [3] , and those for the Zi or Z 2 x Z 2 symmetry in 
the neutrino sector are denoted by Sm, rather than Gj. Following [1], det(T) = det(Siu) = 1 
is assumed, and the order of Wm ■= SnjT is taken to be p. Most of time we will discuss a 
partial Z 2 symmetry and will drop the subscript i in Sjj and Wu- 

The absolute value of a column of the mixing matrix U can be parametrized by six 
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integers (to, k e , k^p, s e , s M ) to be [1] 




or cos ^ + cos -I 2 - — aj sin ^ 





4 sin ^ sin ^ 
a R cos ^ + cos ^ — ai sin ^ 

4 sin sin %^ 
a# cos 4r + cos ^ — a j sin 4p 



[tow] 



(1) 



4 sin sin *f- 



where i labels the column. The integer parameters (to, k e , k^) are defined via the Z m genera- 
tor T = diag(o;^ e , Um , ui^ ), where cu m := exp(27ri/m). The angles are <p e := 2irk e /m, <p^ : = 
2nk fM /m, <p T := 27rk T /m, eM := e -0 M , MT := M -0r, and re := T -0 e . The parameter 
fc r = —(k e + kfj) is determined by det(T) = 1. The quantities or and ai are the real and 
imaginary parts of a := Ti(Wu) = cjp e + oo^ + us s p T , where w* 6 , u^^u^ are the eigenvalues 
of Wjj = S V T. Since det(Wu) = 1, it follows that s e + + s T = 0, so (1) is parametrized 
by six integers (to, k e , k^,p, s e , s M ). The range of these integers are m > 3, m > k e , k^ > 0, 
and p > s e , > 0, but they cannot be chosen independently. 

One indication that the six integers are not arbitrary is the following. The right-hand 
side of (1) is not guaranteed to be non-negative but the left-hand side must be. For example, 
with (to,p) = (3,5), there is at least one negative entry in (1) if we take (k e , k^, s e , s M ) to 



be either (0,1,0,0), (0,1,1,1), (0,1,1,2), (0,1,1,3), (0,1,2,2), or (0,1,2,4), so these choices of 



parameters are forbidden. 

The reason for this failure is that a group structure cannot emerge from these choices, as 
we shall see in Sec. IVB. 

Many known examples in the literature relies on to = 3 and T = diag(l, u, oo 2 ). In fact, 
this T is unique up to permutation for m = 3 if it is to be non-degenerate. For that T, (1) 
is reduced to 



so the equivalent mixing vectors are completely controlled by a = Tr(Wu). We call two 
mixing vectors equivalent if they differ only by a permutation of their components. 




3(^ + 1), 

l(-a R + 2 + >/3aj), 


-(-or + 2 - v^a/), [mv3] 





(2) 
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Table 1 shows some known cases taken from [1-6] with this T. Here the parameters 
(to, k e , kfi) = (3, 0, 1) are fixed. The first column in the Table gives the GAP-designation of 
the group in its Small Group library [7], the second column the common name of the group, 
the third column the mixing vector \v\ — (\U e i\, \ U^i\, \ U T i\), or one of its permutations. The 
last two columns give (p, s e , s M ) and the corresponding (or, aj) parameters that give rise to 
this \v\. (2) can be used to check that the listed \v\ indeed corresponds to the listed (or, ai). 



SmallGroup 


Group 


\v\ 


(P, s e, S M ) 


(a>R,ai) 


[12,3] 


A A 


[.577, .577, .577] 


(3,0,1) 


(0,0) 


[96, 64] 


A(96) 


[.577, .789, .211] 


(8,1,2) 


(0,1) 


[24, 12] 


S 4 


[.816, .408, .408] 


(4,0,1) 


(1,0) 


[150,5] 


A(150) 


[.480, .812, .332] 


(10,1,3) 


(-.309, .951) 






[.777, .607,. 170], 


(10,1,2) 


(.809, .588) 


[162, 14] 




[.525, .279, .804] 


(18,7,13) 


(-.174, -.985) 


[168,42] 


PSL{2, Z 7 ) 


[.408, .894,. 187] 


(7,1,2) 


(-.500,1.32) 


[294, 7] 


A(294) 


[.354, .814, .460] 


(14,1,5) 


(-.623, .782) 


[384, 568] 


A(384) 


[.312, .810, .497) 


(16,1,6) 


(-.707, .707) 


[486,61] 


A(486) 


[.525, .279, .804] 


(18,7,13) 


(-.174, -.985) 


[600, 179] 


A(600) 


[.480, .812, .332] 


(10,1,3) 


(-.309, .951) 



Table 1. Some mixing vectors obtained from (to, k e .k^) = (3, 0, 1) 



III. GROUP STRUCTURE 

Suppose Q is a finite group that contains the elements T, Su, Wjj = SjjT, of orders to, 2,p 
respectively. Then its order \Q\ must be an integer multiple of the least common multiple 
of 2, to, and p. Moreover, the product of an order-2 element Su and an order-m element T 
must produce an order-p element Wu- Whether this can be realized or not depends on Q, 
and is determined by the structure constant of that group. 
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A. Structure Constant 



Let Ci represent the ith conjugacy class of Q, and Cfij the set of all products xy with 
x G Ci and y G Cj. It is not hard to see that Cfij must cover all classes Ck an integer number 
of times so that we can write 



The non-negative integers 7^ are known as structure constants of the finite group Q. It is 
symmetric in i and j because xCjX~ x = Cj, hence QCj = Cfii. They can be calculated from 
the characters \i °f class Ci and irreducible representation ji by the formula 



where is the size of class Cj and is the dimension of the irreducible representation fi. 
Hence {Ckl'jijk is completely symmetric in i,j,k. 

Let c denote the list of classes whose elements have order o. Instead of writing the list as 
{Ci,Cj, ■ ■ ■}, we simply write it as c Q = ■ ■ •}. Given a choice of (m,p), we need to find 
the groups Q that are 'structurally sound', in the sense that there is at least one 7^ > 
with i G C2, j G Cm, and k G c p . A program in GAP language to search for such groups is 
presented in Appendix A. 

B. Eigenvalues 

Given a group Q that is structurally sound for the parameters (m,p), the next task is to 
find out the allowed choice of (k e , k^, s e , s M ) for this group. Since these parameters are related 
to the eigenvalues of T and W v , and all elements within a class have the same eigenvalues, 
they can be obtained by searching the eigenvalues of classes Cj G c m and Ck G c p for which 
lijk > 0. There could be different sets of eigenvalues for different 3-dimensional irreducible 
representations IR3 and all of them should be found. For c p , it is actually only the trace that 
is needed, not the individual eigenvalues. For a 3-dimensional irreducible representation /i, 
a is simply the character Xj- The output of the program in Appendix A contains these 
allowed eigenvalues. 




(3) 



k 




[gabcchi] 



(4) 
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IV. EXAMPLES 



Let us illustrate the discussion of Sec. Ill with some examples. 
A. (m, k e , kfj) = (3, 0, 1), p = 10 

For (m,p) = (3, 10), the least common multiple of (2, 3, 10) is 30, so we need to consider 
only groups Q whose order is a multiple of 30. With det(T) = det(S[/) = det(Wu) — 1, Q 
can be confined to finite subgroups of SU(3), and we will limit ourselves to those with an 
order < 512. From the list of such subgroups [8], only three groups have the right order. 
They are the Small Groups [60,5] =A 5 , [150,5]=A(150), and [300,43]=A(300). 

Table 2 shows these groups, as well as the allowed parameters (s e , s M ) and a = a r + iai 
when the group is structurally sound. This table is taken from the output of the programs 
'gpsearch' and 'gpsu3' in Appendix A. 

The last four columns in the Table display the order-o classes c G , for o — 2, o = m = 
3, o — p — 10, and the list IR3 of 3-dimensional irreducible representations of the group. 
We see that [60,5] is not structurally sound because it lacks order-10 elements. The group 
[300, 43] does have plenty of order-10 elements, but it fails because there is no non-vanishing 
structure constants between classes of orders 2, 3, and 10. That leaves only the group [150,5]. 

A(150) = [150, 5] has one class each of order 2 and 3, four classes of order 10, and eight 
3-dimensional irreducible representations. However, only the combinations in column 4 have 
non-zero structure constants and unit determinants. These combinations are consistent with 
those listed in Table 6 of [6], though they are obtained in completely different ways. The 
allowed parameters given in columns 2 and 3 form two complex-conjugated pairs, with both 
(s e ,s^) = (1,2) and (s e , s M ) = (1,3) appeared in Table 1. According to (2), exchanging a 
with a* only permutes |f 2] and | ^3 1 , so the mixing vectors of (s e , s M ) = (3,8) and (1,2) are 
equivalent, and those of (4,7) and (1,3) are equivalent. 

Any pair (s e , s M ) in the Table can be replaced by an equivalent pair (s a , sp) in which a, (3 
are chosen from e, //, r, because they all give rise to the same a. For example, take the pair 
(s e ,Sp) = (1,2). Since s e + s M + s T = mod p, it follows that s T = 7, therefore any pair 
taken from 1,2,7 will give the same a and the same mixing vector. 

In principle, both s e and can vary from to 9, but according to Table 2, no combination 
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(s e , S/j) other than these four (and their equivalent ones as explained in the last paragraph) 
are allowed, at least not for a group of order < 512. 



SmallGroup 




a = (a R ,ai) 


(i,j,k,IR 3 ) 


c 2 


c 3 


ClO 


IR 3 


[60, 5] 








2 


3 




2,3 


[150,5] 


(1,2) 
(1,3) 
(4,7) 
(3,8) 


u 10 = (.809, .588) 
ujf = (-.309, .951) 
ul = (-.309, -.951) 
uf = (.809, -.588) 


(2, 3, 7, 9), (2, 3, 7, 11) 
(2,3,10, 7), (2, 3, 11,11) 
(2, 3, 6, 9), (2, 3, 10, 5) 

[Z, 6, 11, I ) 

(2, 3, 6, 7), (2, 3, 7, 5) 

(2,3,10,11) 
(2, 3, 6, 11), (2, 3, 7, 7) 
(2,3,10,9) 


2 


3 


6,7,10,11 


4-11 


[300,43] 








3 


2,6 


7-10, 13-18 
22-29, 31-36 


4-36 



Table 2. Allowed groups and allowed parameters for (m,p) = (3, 10) 



B. (m,k e ,k^) = (3,0, 1), p = 5 

Since the least common multiple of (2, 3, 5) is also 30, the groups to be considered here are 
identical to those in Sec. IVA, namely, [60,5], [150,5], and [300,43]. Table 3 is the equivalent 
of Table 2, now for p = 5. We see that [60,5]=A 5 is the only structurally sound group for 
(m,p) = (3,5). 

In principle, s e and s M can range from to 4, but we saw in Sec. II that the combinations 
(s e , s M ) = (0, 0), (1, 2), (1, 3), (2, 2), (2, 4) cannot be allowed because they give rise to a neg- 
ative |t/ Q j| 2 . Let us see why this is so from the group theoretical point of view. According 
to Table 3, the allowed combinations of {s e -,s^) are (1,4) and (2,3), both leading to s T = 0, 
so the allowed pairs should be chosen either from 1,4,0, or 2,3,0, but none of the forbidden 
pairs in Sec. II are one of these, which is why none are allowed. 
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v t > i o 1 1 ( T rni i t\ 
kJllldll \jl (J Liu 


f<? <? "1 


a — \ a Rj a I ) 


J, h i 1SX 3) 








TRo 
irt3 


[60,5] 


(1,4) 
(2 3^ 


-ulil +uj 5 ) = (1.618,0) 
—u)J\ 4- — (—0 61 5 0) 


(2, 3, 5, 2), (2, 3, 4, 3) 
(2 3 4 2W2 3 5 3"l 


2 


3 


4,5 


2,3 


[150,5] 








2 


3 


4,5,8 
9,12,13 


4-11 


[300,43] 








3 


2,6 


4,5,11 
12,19-21,30 


4-36 



Table 3. Allowed groups and allowed parameters for (m,p) = (3,5) 



C. (m,p) = (4,7) 

The programs 'gpsearch' and 'gpsu3' in Appendix A can be used on any (m,p), not just 
m — 3. To illustrate that, consider the example (m,p) = (4,7). In this case, the least 
common multiple of (2,m,p) is 28. Among finite 577(3) subgroups of order < 512, only the 
groups [84,11], [168,42], [336,57] have orders which are integer multiples of 28. Of these three 
groups, only the group [168,42] = £(168) = PSX(3,2) = PSL(2,7) is structurally sound, 
as seen in Table 4. The allowed parameters (s e , s M ) = (1, 2) with s T = 4, and (s e , s M ) = (3, 5) 
with s T = 6, are complex conjugate of each other. The values of (k e , k^, k T ) not listed in 
the Table are both (0,1,3). These assignments are those in known examples [1, 5], what we 
learn from Table 3 in addition is that there are no other combinations of (k e , k^, s e , s M ) that 
is structurally sound for (m,p) = (4, 7), at least not for groups of order < 512. 



SmallGroup 




a = (or, a/) 


(i,j,k,IR 3 ) 


c 2 


c 3 


c 5 


IRs 


[84,11] 








3 




4,9 


4-12 


[168,42] 


(1,2) 




(2, 4, 5, 3), (2, 4, 6, 2) 


2 


4 


5,6 


2,3 






= (-0.5,1.323) 














(3,5) 


U%(1 + OJj + UJj) 


(2, 4, 5, 2), (2, 4, 6, 3) 














= (-0.5,-1.323) 












[336, 57] 








5 


3,7,9,10 


4,20 


4-40 



Table 4. Allowed groups and allowed parameters for (m,p) = (4, 7) 
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D. [600, 179] = A(600) 



We can also use this technique to find all the allowed parameters (m, k e , k^,p, S2, s M ) for 
a given group. As an illustration let us consider A(600) with (m,k e ,k p ) = (3,0,1). The 
allowed values of {p, s e , s M ) and a, together with the equivalent mixing vectors |i>| calculated 
from (2), are listed in Table 5. 

Since the classes of A(600) have orders 1-5, 10, 20, we will consider p = 3,4, 5, 10, 20. 



p 


( Q \ 


a = 


= (a R ,ai) 


\v\ 


Q 
O 


(0, 1 ) 





= (0,0) 


(.577, .577, 


.577) 


A 

4 


in 1 \ 


1 


= (1,0) 


(.816, .408, 


.408) 






-1 


= (-i,o) 


(.000, .707, 


.707) 


5 
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(0, 5) 


-1 


= (-i,o) 


(.000, .707, 


.707) 




(1,2) 




= (.809, .588) 


(.777, .607, 


.170) 




(3, 6) 


M = 


(-.309, .951) 


(.480, .812, 


.332) 




(3, 8) 


-"i = 


(.809, -.588) 


(.777, .170, 


.607) 




(4, 9) 


-u 5 = (-.309, -.951) 


(.480, .332, 


.812) 


20 


(0,5) 


1 


= (i,o) 


(.816, .408, 


.408) 




(1,8) 


ul = ( 


-.809, .588) 


(.252, .800, 


.546) 




(2,4) 


-"1 = 


= (.809, .588) 


(.777, .607, 


.170) 




(2,6) 


M = 


(-.309, .951) 


(.480, .812, 


.332) 




(3,4) 


u 5 = 


(.309, .951) 


(.661, .746, 


.086) 




(4, 13) 


u 5 = 


(.309, .951) 


(.661, .746, 


.086) 




(6,16) 


-"1 = 


(.809, -.588) 


(.777,. 170, 


.607) 




(7, 16) 


ui = ( 


.309, -.951) 


(.661, .086, 


.746) 




(8,18) 


-W5 = ( 


-.309, -.951) 


(.480, .332, 


.812) 




(9,12) 


^(- 


809, -.588) 


(.252, .546, 


.800) 



Table 5. Allowed (p, s e ,s^,a 



and mixing vector \v | for A(600) with (m, k e , k^) = (0, 1, 2) 



A neutrino mixing matrix with the third column given by the mixing vector (.170, .607, 
.177) and the second column given by the mixing vector (.577, .577, .577), both appearing 
in Table 5, is in good agreement with the measured neutrino mixing parameters [2-4, 6]. 
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Appendix A: GAP Programs 

1. Usage 

Two programs written in the GAP language [7] are given in this Appendix. 

The main program is the function gpsearch(m,p,g,LS). It tests whether the structure 
constant 7^ (see (3)) for the group G =SmallGroup(g) is nonzero for some i G c 2 , j G c m , 
and k E c p (see Sec. IIIA). If not, it returns 'structurally forbidden'. If so, it gives different 
printouts depending on whether the short form (LS=1) or the long form (LS=2) is chosen. 

If LS=1, all the allowed [k e , k^, s e , s M , a] (see Sec. II) are printed out. If LS=2, then the 
printout includes the lists c 2 , c m , c p and the list IR 3 of 3-dimensional irreducible representa- 
tions of G. In addition, also a list of k, n], [k e , k T ], [s e , s M , s T ], a] for which 7^ > 0. 
The allowed parameters (k e , k^, k T , s e , s M , s T , a) depend on the irreducible representation n G 
IR3. 

The program stops when one of the following happens. Either the order of G is not a 
multiple of the least common multiple of (2,m,p), or there is no 3-dimensional irreducible 
representation in G, or one of c 2 ,c m , c p is empty. 

The second function gpsu3(m,p) runs gpsearch(m,p,g,l) through all g which is a finite 
SU(3) subgroup of even order less than 512. 

2. gpsearch 

a. Explanations 

• lines 04-10 defines a sub-function which converts a list L of cyclotomic numbers 
E(n) fe := into a list of the exponents k. 

• the lists of C2, c m , c p , and IR3 are computed in lines 24-26 and line 20. 

• the lists dp2, dm2, dp2 in line 35 are sublists of c 2 ,c m ,c p with determinant=+l. 

• the do loop between lines 28 and 59 runs through all n GlR 3 . 

• the three do loops between lines 40 and 50 run through all i G dp2, j G dm2, and k G 
dp2. 
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• 7i?fc is calculated in line 41. 

• the eigenvalues of T are calculated in line 43, and converted to the index form 
[k e , kfj,, k T ] in line 47. 

• the eigenvalues of Wu are calculated in line 48, and converted to the index form 
[s e , Sp, s T ] in line 52. 

• a —Tr(Wu) is calculated in line 50. 

• the rest of the lines prepares for the outputs, depending on whether LS=1 or LS=2. 
b. The Program 



01 gpsearch:=function(m,p,g,LS) 

02 local convert,out,class,n,G,dir,ir3,order,p2,pm,pp,iir3,d2,dm,dp,dp2,dpm,dpp,tbl,ip2, 

03 ipm,ipp,evm,evp,a,kevm,sevp,outl,out2,out3; 
#- 

04 convert:=function(L,n) # L is a list, with L[i]=E(n)Ak[i]. Outputs k[i] 

05 local i,j,k,nl; 

06 nl:=Number(L); k:=List([l..nl],x->0); 

07 for i in [l..nl] do for j in [0..n-l]do 

08 if E(n)Aj=L[i] then k[i]:=j; break; fi; 

09 od; od; return(k); 

10 end;; 
#- 

11 if LSol and LS<>2 then 

12 return("the LS argument has to be either 1 (short) or 2 (long)"); fi; 

13 out:=[ ]; 

14 if not IsInt(g[l]/LcmInt(2,LcmInt(m,p))) then return("group has wrong order"); fi; 

15 # stops the program if lcm(2,m,p) does not divide the order of the group 

16 G:=SmallGroup(g); # G is the group to be searched 

17 class:=ConjugacyClasses(G); # a list of conjugacy classes of G 

18 dir:=List([l.. Number (class)], i— >Trace(Identity(G)AlrrcduciblcRcprcsentations(G)[i])); 
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19 # list of dimensions of irreducible representations of G 

20 ir3:=Positions(dir,3); # a list of 3-dim irreducible representations of G 

21 if ir3=[ ] then return("this group has no 3-dim irred rep"); fi; 

22 # stops the program if there are no 3-dimensional irr rep 



23 order:=List(class, x— 7-Order(Representative(x))); # list of class orders 

24 p2:=Positions(order,2); # a list of order-2 classes 

25 pm:=Positions(order,m); # a list of order-m classes 

26 pp:=Positions(order,p); # a list of order-p classes 

27 if p2=[ ] or pm=[ ] or pp=[ ] then return("this group contains no such orders"); fi; 

28 for iir3 in ir3 do 



29 d2 : =Positions (List (p2 ,x— ^Determinant (Representative ( 

30 class [x] ) AlrreducibleRepresentations(G) [iir3] )) ,1) ; 

3 1 dm: =Positions (List (pm- ^Determinant (Representative ( 

32 class[x])AlrreducibleRepresentations(G)[iir3])),l); 

33 dp: =Positions (List (pp ,x— ^Determinant ( 

34 Representative(class[x])AlrreducibleRepresentations(G)[iir3])),l); 

35 dp2:=List(d2,i— )-p2[i]); dpm:=List(dm,i— >pm[i]); dpp:=List(dp,i— >-pp[i]); 

36 # for each 3-dimensional irreducible representation iir3, 

37 # dp2, dpm, dpp are sublists of lists p2, pm, pp with determinants 

38 if d2=[ ] or dm=[ ] or dp=[ ] then continue; fi; 

39 tbl:=CharacterTable(G); # character table of G used to calculate structure constant 

40 for ip2 in dp2 do for ipm in dpm do for ipp in dpp do 

41 if ClassMultiplicationCoemcient(tbl, ip2,ipm,ipp)=0 then continue; 

42 # select out non-zero structure constants gamma[ip2,ipm,ipp] 

43 evm:=Eigenvalues(CF(m),Representative( 

44 class[ipm])AlrreducibleRepresentations(G)[iir3]); 

45 # eigenvalues of T in that irreducible rep 

46 if Number (evm)<3 then continue; fi; # proceed only if T is nondegenerate 

47 kevm:=convert(evm, m); # convert eigenvalues of T into exponents [ke, km, kt] 

48 evp : =Eigenvalues (CF (p) ,Representative ( 

49 classfipp]) AlrreducibleRepresentations(G) [iir3] ) ; 

50 a:=Trace(Representative(class[ipp])AlrreducibleRepresentations(G)[iir3]); 
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51 # distinct eigenvalues and trace of W_U 

52 sevp:=convert(evp,p); # convert eigenvalues of W_U into exponents [se, sm, st] 

52 outl:=Concatenation([kevm],[sevp],[a]); # list [kevm, sevp, a] 

53 out2:=Concatenation([[ip2,ipm,ipp,iir3]],outl); # list [[ip2,ipm,ipp,iir3],[kevm,sevp,a]] 

54 out3:=[kevm[l],kevm[2],sevp[l],sevp[2],a]; # list [ke,km,se,sm,a] 

55 if LS=2 then out:=Concatenation(out, [out2]); 

56 # string lists out2 together into a big list out 

57 elif LS=1 then out:=Concatenation(out, [out3]); fi; 

58 # string lists out3 together into a big list out 

59 od;od;od;od; # close iir2, ip2, ipm, ipp loops 

60 if out=[ ] then return(" structurally forbidden\n"); fi; 

61 if LS=2 then Print("\nc2=",p2, " cm=",pm, " cp=",pp," IR3=",ir3,"\n"); fi; 

62 if LS=2 then Print (" [ip2,ipm,ipp,iir3], [ke,km,kt], [se,sm,st],a\n"); fi; 

63 return(out); end;; 



3. gpsu3 

01 gpsu3:=function(m,p) # search SU(3) finite subgroups of order <512 

02 local su3fil, g; 

03 su3fil:=[ [ 12, 3 ], [ 24, 12 ], [ 48, 3 ], [ 54, 8 ], [ 60, 5 ], [ 84, 11 ], [ 96, 64 ], 



04 [ 108, 15 ], [ 108, 22 ], [ 150, 5 ], [ 156, 14 ], [ 162, 14 ], [ 168, 42 ], 

05 [ 192, 3 ], [ 216, 88 ], [ 216, 95 ], [ 228, 11 ], [ 294, 7 ], [ 300, 43 ], 

06 [ 324, 50 ], [ 336, 57 ], [ 372, 11 ], [ 384, 568 ], [ 432, 103 ], 

07 [ 444, 14 ], [ 486, 61 ] ];; # even-order finite subgroups of SU(3) 



08 for g in su3fil do if not IsInt(g[l]/LcmInt(2,LcmInt(m,p))) then continue; fi; 

09 Print(gpsearch(m,p,g,l),"\n\n"); od; end;; 
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